Generation of Werner states via collective decay of coherently driven atoms* 
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We show deterministic generation of Werner states as a steady state of the collective decay dy- 
namics of a pair of neutral atom coupled to a leaky cavity and strong coherent drive. We also show 
how the scheme can be extended to generate 2A'^-particle analogue of the bipartite Werner states. 

PACS numbers: 



Entanglement, one of the most striking features of 
quantum physics, is a great resource for the modern fields 
of quantum information and quantum computing 
Entanglement makes possible novel quantum communi- 
cation protocols such as quantum teleportation Q and 
quantum dense coding and is at the heart of quantum 
cryptography 0|. 

Among several bipartite entangled states, Werner 
states |3 provide the simplest example of mixed states 
possessing entanglement. Werner states, in fact, signify 
a wide class of entangled states with a varying degree of 
correlations depending on a single parameter. Assuming 
a general mixed state density matrix pyi for a bipartite 
two-state system, the entanglement fidelity can be ex- 
pressed through F — {^s \ Pm I^'s), where \^s) is a singlet 
state — of the Bell type — of the bipartite system. De- 
pending on this fidelity, the Werner state demonstrates 
classical or Quantum nature. It has been shown by Ben- 
nett et al. jal that for F > 1/2 the mixed state can be 
purified to obtain a pure singlet state, thus can be taken 
to possess non-local character 0. However, for F < 1/2 
it is completely classical and can be generated by using 
an initially unentangled particles. Nevertheless, it has 
been reahzed that only when F > (2 + 3\/2) /8 « 0.78 
the Werner states violate the Clauser-Horne-Shimony- 
Holt inequality Q and demonstrate quantum correla- 
tions. Experimentally, Werner states have been gener- 
ated so far with only photonic qubits f^. Here, for the 
first time we propose a scheme for generation of Werner 
states with neutral atoms. The proposal can be, fairly 
easily, extended to trapped ions and it is well-within the 
realm of current experiments. The scheme is based on 
collective decay dynamics of two atoms coupled coher- 
ently to a strong drive field. As the Werner states are 
obtained as steady state solutions of the decoherent dy- 
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namics, they are naturally stable and long-lived. The 
model can be extended to a 2iV-particle system to obtain 
a multiparticle analogue of the usual bipartite Werner 
state as we discuss in detail. 

To throw some light on the collective dynamics of a 
system of two-level atoms we first consider a simple case 
of two atoms interacting with a leaky cavity with the de- 
cay rate k. The atom-cavity interaction can be described 
through the Hamiltonian 

M' = m {S^ cos ^ + sin () a) + H. c. , (1) 

where is the atom-cavity coupling Rabi frequency. 

The cavity field has the frequency and a is its an- 
nihilation operator. S~ = \gi) {ei\ are the atomic opera- 
tors. The terms cos^ and sin^ are introduced to incor- 
porate the modal dependence of the cavity field and the 
relative position of the two atoms. The density matrix 
dynamics of the complete atom-cavity system is given by 

Pac = ^■^t^' ^ («^«Pac - Sapacfl^ + PacA^a) • (2) 

Using the rotating wave approximation, and after trac- 
ing out the field degrees of freedom the density matrix 
equations for the atoms can be written as 

p=-T {R+R-p - 2R-pR+ + pR+R-) , (3) 

where R^ = cos^ + 5*2" sin^ and F = ^^/k. We have 
assumed that <^ k (bad-cavity limit) while arriving 
at the above equation. The density matrix equation (O 
closely resembles the one describing spontaneous decay 
of an atom with the rate 2r, except that that atomic 
operators R^ are collective and thus it describes collec- 
tive decay dynamics of a system of atoms. The steady 
state solution of the above matrix equation can be readily 
verified to be 

p = D{R-)-\R+)-\ (4) 

provided that the determinant of the operator matrix 
R~ is non-zero and its inverse exists, with D being the 
appropriate normalization. In case the determinant is 



2 



zero, corresponding to the existence of zero eigenvalues 
R- \\Sq) = 0, the solution of Eq. (0) is given by 



P = |*o) (*c 



(5) 



These general considerations hold, immaterial of the ac- 
tual form of the operators R~ and i?"*" , so long as they 
satisfy Lie algebra similar to the spin operators. 

It can be readily verified that the zero eigenen- 
ergy states of R~ are \ipg) = \gi,g2) and |V'e) = 
(cos^ |gi, 62) — sin^ Id, (?2))- The state \tpg) is a sepa- 
rable state and {i/je) is similar in form to one of the max- 
imally entangled Bell-states. It can be noted, however, 
that starting with the initial state lijjg), the entangled 
state \iIje) could not be reached in the steady state. Thus, 
the initial state needs to be asymmetric i.e., with one of 
the atom in the excited state. For example, with the 
initial state \ei,g2), which can be written as 

1^1,62) = a[cos^ 151,62) - sin^ |ei,g2)] 

+ P[sm^\gi,e2)+cos^\e2,g2)], (6) 

with a = cosi^ and (3 — sin^; the steady state density 
matrix is given by 



p = lap IV-e) (M + |/3P 151,52) (51,521 



(7) 



Thus, the entangled state \4>e) is generated with the 
probability |ap = cos^ ^. Nevertheless maximal incoher- 
ent mixing is observed for ^ = 7r/4, 37r/4, 

Our calculations further suggest that even good cavi- 
ties can be used to generate mixed-state entanglement at 
steady state. Here we use the master equation (0 and 
solve it for the time evolution of the density matrix of the 
combined system of the two atoms and the cavity mode. 
To illustrate, starting with an initial state |ei,52), with 
no photons in the cavity mode, mixed-entanglement of 
the form |71l can be generated at the steady state pro- 
vided K 7^ 0. Thus, the ideas presented so far can also be 
applied to two atoms trapped in a good-cavity. 

Thus, we have seen how collective decay dynamics of 
a system of two atoms can give rise to entanglement. In 
the following discussion we add a resonant coherent drive 
field in order to obtain extra handle over the steady state 
density matrix. We assume a general collective-decay 
mechanism which we have shown how to to obtain in a 
leaky cavity for example. 

The dynamics of a system of two level atoms driven 
coherently by a classical field (Rabi frequency |17|e~"'^) 
on-resonance and experiencing a collective decay at the 
rate 2r is governed bv[lo| 



p = i\n\ [e'*5+ + e-'*5-,p] 
- T{S+S-p - 2S-pS+ + pS+S- 



(8) 



Here, are the collective atomic operators; they can 
be thought of as the operators for a spin N particle for 



a total of 2N particles in the system. These collective 
operators are very similar to the regular spin operators 
except for the absence of h, dependence of their eigenval- 
ues unlike the actual spin systems. 

Introducing a new operator defined as 
R- = S- +i{\n\/T)&''' m [13, the density matrix 
equation ((SJ takes the same form as Eq. JH)). Therefore, 
the steady state solution, which is analogous to Eq.Q), 
can be readily arrived at. In the strong-field limit with 
an initially symmetric state of the atomic spins the 
steady state solution is given by 



p(oo) = 



1 



25 + 1 



s 

E 

n=-S 



S, m) {S, TO I 



(9) 



which is a mixed state with all possible eigenstates 
equally occupied. Moreover, the state \S = 0,m = 0) is 
the zero eigenvalue state giving a pure state solution if 
the initial state overlaps with it. 

In the following, we discuss how the collective decay 
dynamics can be exploited to generate a bipartite Werner 
state of a pre-chosen fidelity. We start with a general ini- 
tial state sm9 |ei,52) -|- cos9 151,62), for a bipartite sys- 
tem. It can be readily rewritten as 



sm ( 



V2 



|*+) 



sm f — cos ( 



V2 



(10) 



in terms of two of the Bell States 



|^±^ ^ leiea) ± I5152) ^ I6152) ± \g1e2) ^^^^ 



V2 



V2 



The Bell states can be readily written in terms of the 
collective spin states — \S, to) — as 

|vl/+)^|l,0),|<l>±)^-i=(|l,l)±|l,-l)). (12) 

Whereas, the state l^*^), which can be written as 
IS* = 0, TO = 0) in the collective spin description, does not 
evolve under Eq. ||HJ). 

Therefore, to study the dynamics of a given initial state 
it is important to decompose it into the states | S, m) 
(with 5 = or 1) or alternatively Bell states as discussed 
above. The density matrix of such an initial state can be 
represented as 



p{t = Q) = F\^-) 1 + (1 - i^) pT . 



(13) 



Thus, for the initial state (|10|l . the fidelity of the maxi- 
mally entangled singlet state is 



F= [1 - sin(26l)] /2 . 



(14) 



To recall, |^^) being the singlet state does not evolve 
under the action of ijSJ), and the triplet state density- 
matrix pt evolves with time; and in the strong drive 



limit we obtain 
1 



/3t(oo) = 



E 

m=-l,0,l 



1, m) (1, m| 



(|^+) (^+1 + |$+) ($+1 + 1$-) ($-|) . (15) 



Thus, the coUective decay dynamics of a strongly 
driven two particle system initially in the state 
sin6' Id, o^} +COS0 Igi, 62) gives the Werner state (see for 
example p): 



jl-F) 



3 

- |$+)($^ 



(16) 



with _F as the pre-chosen probability of the singlet state 
component. 

To quantify the parameter 0, we note from Eq. H14|l 
that e = (1/2) sin~^(l-2i^). To obtain purifiable Werner 
states with fidelity -F > 1/2 it is required that 6 e 
(7r/4 + n7r, 37r/4 + n7r) where n is a non-negative integer. 

The time evolution of the non-singlet component of the 
initial state 1)10(1 (a) is shown in Fig.^to show that it ap- 
proaches the non-singlet component of the Werner state 
in Eq. I|16|) . The von Neumann entropy of the Werner 
state ifT^ is given by 



= FlnF+(l-F) ln(l-F)+ln(l/3) . (17) 



^Werner 



Whereas, the entropy of the general steady state density 
matrix as a function of a finite parameter fl/T governing 
the dynamics given by Eq. of the initial state Eq. H13|) 
can be shown to be 

S{t^oo) ==FlnF+(l-F)[ln(l-i^)+/3(r2/r)j (18) 

where /3(f7/r) = J2^ Pi,iHPi,i) fori = {-l,0,l}^he 
steady-state density matrix elements are given by |l3l | 

pia = (l,l|p|l,l> = |xlV^ 
Po^o^{l,0\p\l,0)=-\x\'{l-\x\')/D 
P-u-i = (1, -1| p |1, -1) - 1 - pi,i - po^o (19) 

where x = iV2r2/r and D = Z\x\^ ~ 2|xP + 1. The 
function /3(r2/r) is plotted in Fig. ^(b). It approaches 
ln(l/3) so that the entropy 118(1 matches with that of the 
Werner state given in Eq. 1(17(1 as fi/F ^ 1. 

It can be noted that the proposed scheme is very gen- 
eral and can be applied to multipartite systems as well. 
We start with a general 2iV-partite density matrix 



p(t^^)^F\^,) {^^\^(\-F)p^. 



(20) 



Here — 15 = 0, m = 0) and for the initial density 
matrix devoid of the 5* = component we assume a very 
general form given by 



N 



PT{t = 0)^Yl E E c.^^L\S,m){S,m'\. (21) 




FIG. 1: (a) Time evolution of the non-stationary compo- 
nent of the initial state Eq. (O, i.e., (*+|p|*+) /(I - F); it 
approaches the steady state value 1/3 for sufficiently strong 
drive Q/F = 5 along with the other components ($^|p|<I>^) 
(not plotted here) as given in Eq. Ill6t . (b) The non-trivial 
term in the entropy equation (1181 — f3{Q/T) — of the steady 
state. It approaches to that of the Werner state — ln(l/3) — as 
n/r 2> 1 shown by the dashed line. 



It needs to be noted at this point that the dynamical 
equations conserve the spin value S. Thus, each spin- 
multiplet (m ~ — 5, — S* + 1, • • • ,S) evolves on its own 
without any coupling to other spin counterparts possible 
for the 27V particle state and approaches the state given 
in Eq. (O at the steady state in the strong field limit. 
Thus consolidating all the spin components we obtain 
the steady state density matrix as 

p{t = cx,)=F\-i's) (^'.| + (1-F)x 

N S 

- \S,m){S,m\, (22) 



E 



25 



where a^^'^ = J2m=-s |c*m,mP- The above state can be 
recognized as a 27V-particle analogue of the Werner state 
representing an entangled mixed state parametrized 
by the fidelity of the maximally entangled state. 
Now we use the angular momentum algebra, such 
that for four particles the individual spin vectors 
are added in a manner depicted by the notation 
|5i, ^2, [5'], 53, S4, [S"];S,m = mi + m2 + + m^). 
Thus, the angular momenta of the first two and last 
two particles are added together to obtain S' and S" 
and then combined to obtain total angular momentum 
quantum numbers S and m. In general, there are 
several other ways of adding the given number of 
individual angular momentum vectors, but they are not 
independent and can be transformed into one another 
by simple linear transformations. However, consistently 
using only one such addition scheme is important in 
one calculation. Using the above addition scheme the 
various 4-particle spin states can be written out in terms 
of the bare atomic states as: 



S—l m— — S m' — — S 
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|2,0) = ( |ei,e2,53,54) + |ei, 32, 63, 34) + 131,62,63,54) 
|0,0)i EE (2 161,62,53,54) - [61,52,63,54) - 151,62,63,54) 
10,0)2 = (|ei,52,63,54) - 161,52,53,64) - 151,62,63,54) H 

|1,0)^ EE (|ei,62,53,54) - |5l,52,63,64))/\/2 

11,0)2,3 = (|ei,52,63,54) T 161,52,53,64) ± 151,62,63,54) - 



+ 161,52,53,64) + 151,62,53,64) + 151,52,63,64)) /Ve 

- |ei,52,53,e4) - 151,62,53,64) + 2 I51, 52, 63, 64)) /(2V3) 

- 151,62, 53, 64)) /2 

- 151,62,53, 64)) /2 (23) 



It can be noted that the spin states |0, 0) and |1,0) 
possess multiple representations due to the degeneracies 
present in the addition of several angular momenta to ob- 
tain same resultant. Thus, the state \S — 0,m — 0) can 
be obtained by two different angular momentum addition 
paths such as |0, 0)^ = |l/2, 1/2, [1], 1/2, 1/2, [1], 0, 0), 
and 10,0)2 = |l/2, 1/2, [0], 1/2, 1/2, [0], 0,0). Sim- 
ilarly, |1, 0)1 = |l/2, 1/2, [1], 1/2, 1/2, [1], 1, 0), 
11,0)2 - |l/2, 1/2, [1], 1/2, 1/2, [0], 1,0), and 
|1, 0)3 = 11/2, 1/2, [0], 1/2, 1/2, [1], 1,0). The set in 
Eq. (|23|l . with trivial additions of the states not shown, 
is orthogonal and can be used to uniquely represent any 
given atomic state. 

We consider a 4-particle initial state, |V'4(0)) = 
sin 6* 1 61, 62, 53,54) + cos 9 |5i,52, 63, 64), and show gener- 
ation of a generalized Werner state at steady state. This 
state can be rewritten as 



entangled states of a general 2N particle system which 
would be essentially multiparticle analogue of the usual 
bipartite Werner states. To our knowledge, this is the 
first proposal showing generation of Werner states with 
atomic qubits. We further expect to see similar results in 
other contexts such as collective dephasing in two qubits 
in quantum-dot systems 
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|V'4(0)) =sin( 



-i=(V2|0,0)i + |2,0)-fV3|l,0)i 



cos 9 



^(V2|0,0)i + |2,0)-x/3|l,0)i 



(24) 



Thus the initial entanglement fidelity is given by = 
(1/3)[1 -I- sm(29)]. which is the same as the maximum 
entanglement fidelity possible for the generalized Werner 
state for four particles. This initial 4-particle state, at 
the steady state and in the strong drive field limit, would 
take the from of Eq. ^ with TV = 2, 



^ 3 [1 ~ sin(2g)] ^(2) ^ 1 [1 + sin(2g)] 

2 [2 - sin(26')] ' 



2 [2 - sin(26l)] ' 
and F = [1 + sin(26i)]/3 



(25) 



generating a 4-particle Werner state with \'^s) = |0, 0)-^. 

In conclusion, we have shown how collective decay 
dynamics of the system of two level atoms interacting 
with a leaky cavity develops giving rise to entanglement 
at steady state. With addition of an extra control, in 
the form of a strong coherent drive, the dynamics shows 
further interesting features allowing generation of bipar- 
tite mixed entangled states, namely Werner states. The 
technique can be easily extended to generation of mixed 
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